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The dynamics of entanglement has recently been realized as a useful probe in studying ergodicity 
and its breakdown in quantum many-body systems. In this paper, we study theoretically the growth 
of entanglement in quantum many-body systems and propose a method to measure it experimentally. 
We show that entanglement growth is related to the spreading of local operators in real space. We 
present a simple toy model for ergodic systems in which linear spreading of operators results in 
a universal, linear in time growth of entanglement for initial product states, in contrast with the 
logarithmic growth of entanglement in many-body localized (MBL) systems. Furthermore, we show 
that entanglement growth is directly related to the decay of the Loschmidt echo in a composite 
system comprised of several copies of the original system, in which connections are controlled by a 
quantum switch (two-level system). By measuring only the switch’s dynamics, the growth of the 
Renyi entropies can be extracted. Our work provides a way of understanding entanglement dynamics 
in many-body systems, and to directly measure its growth in time via a single local measurement. 


INTRODUCTION AND RESULTS 

In the past decade, quantum entanglement has 
emerged as an indispensable tool for characterizing and 
classifying ground states of many-body systems, for ex¬ 
ample in the field of topological order m- 

Recently, it was realized that entanglement dynam¬ 
ics also exhibits universality, providing a useful probe 
of ergodicity and its breakdown in many-body systems. 
In particular, the growth of entanglement following a 
quantum quench can be used to distinguish between lo¬ 
calized and ergodic phases. For many-body localized 
(MBL) systems [3, 5], emergent local integrability 0- 
[5 implies that entanglement grows logarithmically in 
time, S(t) ~ logf, following a quantum quench from ini¬ 
tial product states [9] [10]. For generic quantum ergodic 
systems, entanglement appears to grow universally lin¬ 
early [TTj: S(t) ~ t. 

Previous works predicted linear growth S(t) ~ t 
in integrable systems, including (1 + 1)-dimensional 
CFTs |T2j 1 1 ■I| ■ while other works studied entanglement 
growth in non-integrable, higher dimensional CFTs via 
holographic calculations ® na. However, an under¬ 
standing of the entanglement dynamics in generic quan¬ 
tum chaotic/ergodic systems is lacking. One puzzling 
fact pointed out in Ref. m is that energy transport in 
a quantum chaotic or ergodic many-body system under 
a time-independent Hamiltonian is diffusive (~ y/t), yet 
entanglement growth is linear (~ t) - clearly, the mech¬ 
anisms of particle and quantum information transport 
are different, and so it is important to understand the 
dynamics of the latter. 

The purpose of this paper is to provide 1) a theoretical 
description, and hence, a physical picture of the growth 
of entanglement entropy (EE) in a many-body system, 
and 2) a proposal to experimentally measure it. To be 
precise, we consider the dynamics of the ro-th Renyi EE 


S n (t) of pure states | ip) that are initially random product 
states, evolving unitarily under a local, potentially time- 
dependent non-integrable many-body Hamiltonian: 

H = Y j Hx- ( 1 ) 

The Hamiltonian is defined on a lattice A of sites (labeled 
by i ) in d spatial dimensions so that the local Hilbert 
space Hi is bounded: dim('Hj) = k < oo, and X is a local 
region in A. The n-th Renyi EE of A , for a bipartition 
of the system into two subregions A and B , is given by 

S n (t) = —— log Tr(p^(t)), (2) 

1 — n 

with pa{1) = Tr s (t/ t |'0)('i/>|t//) being the reduced density 
matrix of subsystem A , and Ut the unitary time evolution 
operator generated by H according to the Schrodinger 
equation. The state, being initially a random product 
state, has S„(0) = 0, but S n (t) > 0 for t > 0 and grows 
in general. Also, S n > S m for any m > n\ in particular, 
S 2 is a lower bound for S, the von Neumann entropy. 

A summary of our results are as follows. We show that 
the growth of entanglement as measured by S? (t) is di¬ 
rectly related to the measurement of basis operators of 
subregion A. Under time evolution by a local Hamilto¬ 
nian, a basis operator physically spreads in real space, 
and as it grows, the value of its measurement typically 
decreases, thus leading to entanglement growth. Fur¬ 
thermore, we introduce a simple toy model for ergodic 
many-body systems in which 1) operators spread at some 
maximal velocity v, 2) delocalize completely within this 
light cone r ~ vt, and show that ^(t) grows linearly 
in time with a velocity related to v. We believe such a 
model captures the salient features of the universal linear 
in time growth of entanglement seen by Ref. [11] . 

We also propose a way to experimentally measure S n (t) 
via a local measurement. We introduce a quantum switch 
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FIG. 1. Set-up to measure the growth of the n-th Renyi 
entropy S n (t). Here n = 3 and we have shown it for a d = 
1 chain. There is a quantum switch in the middle of the 
set-up which governs tunneling between different subsystems, 
(a) When the quantum switch is in the state | f), tunneling 
between A; and Bi is allowed while tunneling between At and 
Bi +1 is prohibited, (b) When the quantum switch is in the 
state | J.), the allowed and prohibited tunnelings are swapped. 
A composite state in an n-copy product state of the chains 
and in a superposition of | f) and | 4-} of the quantum switch 
will have two parts evolving in time differently according to 
the quantum switch; measuring a x (t) of the quantum switch 
gives the Loschmidt echo and hence the n-th Renyi entropy. 

(a 2 -level system) that allows tunneling between differ¬ 
ent parts of a replicated system consisting of n disjoint 
copies of the original system, depending on the state of 
the switch m- By preparing the replicated system and 
the quantum switch appropriately, and by subsequently 
measuring a x ( t ) of the quantum switch only, the entropy 
growth S n (t) corresponding to the entanglement entropy 
of the original system can be measured, see Fig. [l] 

We will refine and expound upon these ideas below. 
For clarity of argument, we will focus on the case of a 
spin- 1/2 system, so dim (Hi) = 2 , but our results are 
general and hold for systems with bounded local Hilbert 
spaces of other dimensions m We also choose subregion 
A to be a ball of radius r A with volume jdf A = Na in 
d dimensions (7 d = r(I/ 2 +i) the volume of a unit d- 
ball), but the analysis can be readily adapted to other 
geometries. 

LOCAL MEASUREMENTS OF BASIS 
OPERATORS GIVES ENTANGLEMENT 

We first show that the second Renyi entropy S 2 (t) is 
directly related to measuring all time-evolved basis oper¬ 
ators 0 A (t) = UlOAUt in the initial state | i/j), such that 
the operator Oa at t = 0 has support strictly in A (i.e. 
where it acts non-trivially on the lattice). That is, for 
a spin-1/2 system of Na sites, S 2 (t) = — logTr^/X^f)), 
with 

Tr A (paW) = (3) 

o A 


where the sum goes over all operators Oa with support 
strictly in A , assumed to be Hermitian and indepen¬ 
dent as defined by the Hilbert-Schmidt inner product: 
^/jTr a{ 0\0' a ) = So At o> A , and which form a basis for 
all Hermitian operators on A. A potential basis set is 
{I-i, erf , of, of }, where I i,a“ are the identity matrix 
and Pauli-matrices acting on site i respectively, and we 
will use this basis for our subsequent analysis. The proof 
of this statement is straightforward - one notes that the 
reduced density matrix p A (t) = Tr^ \ij>(t)}(ip(t)\ is in par¬ 
ticular a Hermitian operator, and so must be a linear 
combination of Oa operators with real coefficients given 
by 2wxWt)|0A]V’( i )} = Then, squar¬ 

ing p A (t) and using the orthonormality of Oa under the 
inner product, one obtains the claimed result. 

PHYSICAL SPREADING OF BASIS 

OPERATORS LEADS TO ENTANGLEMENT 
GROWTH 

We now use the above result to understand why en¬ 
tanglement typically grows in many-body systems. We 
consider the ensemble average of S 2 (t) over initial pure 
product states \ijj) = \cf\(j 2 - ■ - ctjv ), where 0 / is the ran¬ 
dom direction that the ith spin is pointing to on its Bloch 
sphere. 

It is convenient to switch from an average over an en¬ 
semble of initial states which are random product states 
to an average over an ensemble of locally-rotated Hamil¬ 
tonians. A given initial product state can be written as 
|cri<7 2 • • • a N ) = IL Vi\ tit 2 •'' tiv) where V) is the local 
unitary that rotates | fi) into cq). We can rewrite (|3j) as 

Tr A{p\{t)) = £<tn \U^O A VU t \ ttt} 2 

Oa 

= (4) 

Oa 

where V = JX, Vi and U t is the time evolution opera¬ 
tor generated by the locally rotated Hamiltonian H = 
V^HV. This local rotation generates a new Hamilto¬ 
nian H that has the same locality properties as the 
original Hamiltonian H, i.e. H = which im¬ 

plies that both H and H have the same Lieb-Robinson 
(LR) velocities vlr [T8H2D] . bounds which govern the 
speed of information and operator spreading in quantum 
many-body systems. The second equality arises from a 
straightforward statement of ‘basis invariance’ of mea¬ 
surements B23- Thus, taking an ensemble average of 
S 2 (t) over initial product states in ^ is equivalent to 
taking an ensemble average over locally rotated Hamil¬ 
tonians H in ([ 4 ]) with the measurement done in the par¬ 
ticular state | ttt) - 

We are now in a position to gain a physical understand¬ 
ing of how entanglement grows in a many-body system: 
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under time evolution, a basis operator Oa with support 
in A spreads on the lattice to become Oa (i) = UjOAUt 
(note the evolution under H ), which has not only larger 
support but is also a complicated sum of other basis op¬ 
erators, though with conserved ‘weight’, for example, 


uh 


c 77 — n x rr x 4 - r X ' V 
U t — c i a i + C- 


■»,i+1 L 


T X rry 


’i +1 


Z,X,Z Z 

C i-l,i,i+l a i-l C 


J i +1 


= £< 

X,fi 


X U X j 


(5) 


such that the total ‘weight’ Xx/i( c x) 2 = 1 for 
times El- Then, the value of the measurement (ftt 
|0 a(£)| ttt) typically decreases, as the complicated sum 
of operators will have many ‘off-diagonal’ operators such 
as af of +1 etc. that do not contribute, when only the ‘di¬ 
agonal’ operators which are products of erf s do, leading to 
an increase in entanglement entropy as calculated by ([4]) . 
Thus we see that entanglement growth in a many-body 
system is intimately related to the physical spreading of 
basis operators in real space: in some sense, EE increases 
because quantum information is ‘lost’ in the inability of 
the state | ftt) to measure the increasingly complicated 
operator Oxit). 


A TOY MODEL FOR EXPLAINING UNIVERSAL 

LINEAR GROWTH IN ERGODIC SYSTEMS 

Equations © and © are true regardless of the sys¬ 
tem in question, be it ergodic or otherwise. The differ¬ 
ences in entanglement growth between different systems 
are completely captured in how an operator spreads and 
is decomposed in terms of other basis operators, c.f. |5]). 
Here, let us introduce a simple toy model where operators 
completely ‘scramble’ in a linear lightcone, which gives 
a linear growth of entanglement. We believe that such 
a model captures the salient features of universal linear 
entanglement growth in generic |21j ergodic many-body 
systems, as seen by Ref. El¬ 
ite local Hamiltonians with a bounded local Hilbert 
space, the velocity at which basis operators spreads can 
be at most linear, with an upper bound given by v lr- 
In other words, the operator can only spread within the 
LR light cone r ~ VLRt■ The precise distribution of the 
coefficients {c({} in ([ 5 J) depends on the Hamiltonian in 
question. In our toy model, we make a statistical state¬ 
ment about the distribution of coefficients. Let us assume 
that because of ergodicity, 1) a basis operator spreads lin¬ 
early at some fixed velocity v < Vlr (for all H) so that 
an initially local operator at i spreads to become a sum of 
operators contained within a ball of radius (vt) d centered 
at i, and 2) it is effectively ‘scrambled’ in this ball. Pre¬ 
cisely, we assume that an initially localized basis operator 
at i has a decomposition as in ([5]) under time evolution, 
such that {c£} is a unit random vector of 4 7 d 0 *) coef¬ 
ficients, where X, the support of each basis operator in 



FIG. 2. (Color online) (a) The support of a typical basis op¬ 
erator Oa- Here A is a d = 2 square lattice, and the subregion 
A is a circle of radius va demarcated by the blue circle. Black 
sites represent the presence of a non-trivial operator, erf, erf 
or erf, while white sites represent the presence of an identity 
operator Ij. One can see clusters of white sites within the 
circle. A site within a cluster is ~ 0(1) sites away from a 
black site, (b) Under time evolution, each black site spreads 
with velocity v, and so the clusters quickly get filled up (red 
sites) in time vt ~ 0(1), and simultaneously the operator also 
grows in physical size to become a ball of radius va + vt. 


the decomposition is contained entirely within a ball of 
radius (vt) d centered at i. Such a picture is indeed sup¬ 
ported by Ref. [22], and we assume this behavior to be 
generically true for ergodic systems. We note that such a 
scrambling assumption, tied to ergodicity, will manifestly 
not hold for MBL systems: basis operators that have sig¬ 
nificant overlap with the local integrals of motions stay 
localized near site i, invalidating the assumption. 

Next, we estimate the measurement (ftt \Ox{t)\ ftt) 
in Q. In principle, we have to sum over all operators 
Oa, but we can replace all Oa -a Ox,typ., where Oa, t yp . 
is a typical operator in subregion A having support of 
size L = 3/ANa- (Recall that the support of an op¬ 
erator are the sites where Oa acts non-trivially on, for 
example, Oa = erf <8> I 2 ® of ® I 4 0 • • • <3 Iw A has sup¬ 
port on sites 1 U 3 and size 2.) The statement of the 
size being L = 3/ANa follows because in constructing 
some Oa, there is a choice of three non-trivial operators 
{uf, erf, erf } compared with a choice of a single trivial op¬ 
erator Ij for every site, so there are ( N ^/)3 L operators with 
support of size L , a skewed binomial distribution, which 
gives the average size L = 3/ANa- Typicality follows be¬ 
cause the standard deviation of the size of an operator 
~ \JNa -C Na, coming from the same skewed binomial 
distribution. Thus, Oa( t) can be replaced by Oa, t yp .(£) 
because the latter dominate the sum. 

We now have to understand the time-evolution of 
Oa, typ.- Considered as an operator on the subregion A, 
Oa, typ. has Na — L sites on it acts trivially, i.e., with the 
identity operators. We call groups of contiguous IjS ‘clus¬ 
ters’. A typical cluster has size ~ 0(1). Over time, each 
non-trivial single-site operator that makes up the typical 
operator spreads, emitting a spherical ‘wavefront’, and 
so the clusters quickly get filled up after vt ~ 0(1). Af- 
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ter such time, the resulting operator will therefore have 
grown to become a complicated operator within a ball 
of radius (r A + vt) d , see Fig. [ 2 ] Following the previous 
discussion, we assume that the vector of coefficients {c^-} 
denoting the decomposition of a typical time evolved op¬ 
erator into the basis operators of the ball is also a random 
vector but now with dimension 4 T d(rA+vt) _ 

Next, we estimate the measurement of a typical oper¬ 
ator. This is given by 

07 d(r A +vt) d 

(nt|o«„,(t)im) 2 = SFIW . (6) 


which can be understood as follows. In the decomposi¬ 
tion e>A,typ.(t) = U}0 A ,typ.U t = X)oC 0 (t)C>, the only 
operators O which give a non-zero overlap with | j'ft) 
are the ‘diagonal’ operators I, of, of crj, ■ ■ •. There are 
27 d(r A +vt) d operators contained in a ball of ra¬ 

dius (r A + vt) d . The coefficient co that accompanies 
such an operator is a random variable (which follows 
from our scrambling assumption) and has typical magni¬ 
tude l/y / dim(Ball)(t) = 1/V4A d ( rA+vt ') d . Furthermore, 
since the coefficients co ( t) are independent random vari¬ 
ables for different operators O, the cross-terms in the 

LHS of (6), (Eo c o( 0) ) vanish, reducing the sum to 

Eo co(t) 2 = , the RHS. 

Finally, we arrive at the expression for S 2 (t) and obtain 
the universal linear growth. Plugging the above result 
into Q and simplifying, we get 

Tr A (p A {t)) = 2-^ rA+vt ^- rd ^. ( 7 ) 


This calculation is valid until vt < (2 1 ^ d — 1 )r A , at which 
point Tr ’ A (p A (t)) ~ 1/2 Na . Then, the error from ne¬ 
glecting the contribution of the global identity term 
becomes non-negligible, since the identity operator never 
spreads under time evolution and will always contribute 
a factor of 1/2 Na . Thus, the EE saturates to the maxi¬ 
mum allowed for region A. The second Renyi EE is thus 

s (t) = / 'Yd(i rA + vt) d — r A ) log 2 for ^ <t < 

2 \ld.r d A log 2 = N a log 2 for t> 

( 8 ) 

where fd = (2 1 / d — 1). For times t < 2r A /v{d— 1), S^it) 
can be approximated by a linear function 

Si(t) ~ [( 7 ddr^^vt] log2. (9) 

In fact, this linear approximation is always valid as long 
as the S‘ 2 (t) calculation is valid, since 2 r A /v(d — 1) > 
fdi" A /v for all d. We see that the expression for the dy¬ 
namics of entanglement, ([9]), has three terms: 1) a geo¬ 
metric factor jddr^ 1 which gives the area of an entan¬ 
glement ‘wavefront’, 2) the speed of the wavefront v, and 
3) a universal linear time dependence t. Thus this has the 
interpretation of entanglement spreading in a ‘tsunami’, 
similar to a view put forth in Ref. [SI- 


PROPOSAL TO EXPERIMENTALLY MEASURE 
SJt) AND RELATION TO LOSCHMIDT ECHO 

Here we propose a way to experimentally measure 
S n (t) via a local measurement. We consider an repli¬ 
cated system consisting of n disjoint copies of the orig¬ 
inal system, each with a similar bipartition Ai and Bi. 
We arrange them in a way similar to Ref. |16| : a star 
geometry such that the boundaries between subsystems 
A, and Bi are placed near each other, see Fig. [I] We 
then introduce a quantum switch (a 2-level system), 
which if in the | f) state, allows tunneling only be¬ 
tween subregions Ai and Bi such that the full Hamil¬ 
tonian on the composite system is H , while if in the | j,) 
state, allows tunneling only between subregions Ai and 
Bi + 1 , such that the full Hamiltonian is H + V. Here 
II = E "=1 Hi = E "=1 Ha, + II n, + H Ai B i, where each 
Hi is a copy of the Hamiltonian, ([!]), acting on the zth- 
copy of the Hilbert space Hi = H Ai 0 H n-, and V the lo¬ 
cal ‘reconnection’ operator, V = Y^ 7 i=i H Ai B i+1 — H Ai B t , 
with the cyclic condition n + 1 = 1, so that H + V = 
ELi H,Ai + Hsi + H Ai s i+ 1 - If we prepare the state of 
the replicated system as \i/j n ) = where \ip) is 

the initial product state of the original system and also 
the quantum switch as a maximally entangled Bell state, 
so that the full state is |’F) = \ip n ) <g> -4=(|t) + |4_)), then 
under unitary time evolution, 

|T(t)) = T ( e~ im \ij> nj [) + e ~^ H+v)t \. (10) 

Measuring a x (t) of the quantum switch, a local measure¬ 
ment, gives 

Tr M(t)) = (r\e i{H+V)t e~ iHt \r) = T(t), (11) 

from which the nth Renyi entropy can be obtained: 
S n (t) = yEt logJ 7 )!). Thus we see that the dynamics 
of entanglement is indeed an observable local quantity if 
we work in an extended system. See also Refs. [251121] on 
other theoretical proposals to measure EE and its growth. 
Additionally, Ref. [25] has experimentally measured EE 
in a Bose-Hubbard system also in a replicated set-up. 
However, the difference of our proposal with their exper¬ 
imental technique is that we only require a single local 
measurement wherelse [251 requires an extensive number 
of measurements (namely, measuring the parity on all 
sites). 

The above claim is based on an alternative reformula¬ 
tion of S n (t), relating it to a Loschmidt echo F(t) [25]- 
|5I] on the replicated system. It is a well-known trick 
[T71 1 2511521 that the trace of the nth power of p A can be 
calculated by replicating n copies of the original Hilbert 
space H = H A 0 Hb, and calculating the swap opera¬ 
tor S that cyclically swaps subregions A t —> A imo( i n +i- 
if Q = p® n is the n replicated density matrix, then 

TrM) = Tr H9n (SlS). 
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Now let us specialize to the case considered before: 
when the initial state on one system is a random prod¬ 
uct state I'i/') = \a 1 a 2 ■ ■ -aw) and p = \ip)(i/j\. If each p 
evolves independently under the Hamiltonian given by 
0 . then fl evolves under the composite Hamiltonian 
H = £” =1 H., Then. 

Tr = (V’l e im Se~ iHt \^ n ). (12) 


However, because the initial state is a product state, 
<^”|S+ = (ip n |, and since &e iHt S = e *( ff +v)*, we end 
up with the claimed result, (11), which is a Losclimidt 
echo J-(t) with ‘perturbation’ V to H. Note that an anal¬ 
ogous statement also holds for time-dependent Hamilto¬ 
nians even though the unitary time evolution operator 
will not have the form e~ lHt . 


CONCLUSION 

We have provided a theoretical description and a phys¬ 
ical picture of the entanglement dynamics in a many- 
body system by relating entanglement growth to the 
physical spreading of basis operators on a lattice. We 
also introduced a simple toy model where an initial ba¬ 
sis operator effectively scrambles in its light cone, thus 
producing a linear entanglement growth, which we be¬ 
lieve captures the salient features of the universal linear 
in time growth of entanglement seen in ergodic many- 
body systems. Furthermore, the entanglement growth 
in this model has an interpretation of an ‘entanglement 
tsunami’: there is a prefactor which depends on the ge¬ 
ometry of the subregion A, which propagates outwards in 
a ‘wavefront’ at speed v. We remark that while entangle¬ 
ment growth in MBL systems has been successfully ex¬ 
plained through a local integrals of motion picture, it is 
an interesting question to explain this dynamics in terms 
of our operator spreading language. 

We have also provided an alternative interpretation of 
the growth of the n-th Renyi entropies as a Loschmidt 
echo in a composite system, subject to a perturbation 
that reconnects different subregions. By using this re¬ 
formulation, we have proposed an experimental way of 
measuring the growth of the Renyi entropies: one can ef¬ 
fect the reconnection by a quantum switch and measure 
the state of the quantum switch to extract the Loschmidt 
echo and hence the EE. This proposal can be used, for 
example, to directly detect MBL phases by measuring 
the logarithmic growth of entanglement. 

Note added. - Recent works [(IT 133] also studied en¬ 
tanglement dynamics in chaotic many-body systems, and 
argued that entanglement growth speed is bounded from 
above by the operator spreading speed. In our operator 
counting language picture, this can be accounted for in 
our toy model by refining the ‘scrambling’ assumption 
which gives rise to the measurement value §: there will 
be a distribution of coefficients that presumably depends 


upon the size of the operators in the decomposition ([5|, 
which we leave for future work to explore. 
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is, 


dr a log k = N A \ogk 


^d{{r a + vt) d - r d A ) log k for 

for t > 

(15) 
S 2 (t) 


where fd, = (2 — 1). For times t < 2r A /v(d — 1) 
is approximated by a linear function 


S 2 {t) ~ ((7 d dr d A l )xvxt) log k 


(16) 


which is always valid as long as the S 2 (t) calculation is 
valid, since 2r A /v(d — 1) > fd.r A /v for all d. We see 
that the universal linear growth of EE is still present for 
systems with different dim('Hj) = k < 00 . 


II. Derivation of results 


Here we perform two things. First, we prove (13), that 


Appendix A: Analysis for systems with bounded 
Hilbert spaces of other dimensions 


Tr A (p\(t)) = -E E (Tr {pO A (t)))\ (17) 
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Here we consider the extension of our analysis and re¬ 
sults for systems with bounded Hilbert spaces of other 
dimensions, dim(7f i ) = k < 00 , so that the full Hilbert 
space is H = < 8 > 


and second, we prove (14), 


Tr A {p\{t)) = k-^ rA+vty 


~r d A ) 


(18) 


repeating the arguments presented in the main text for 
a system with dim ('Hi) = k < 00 . 


I. Results 


Proving (13) 


We simply state the corresponding results here. We al¬ 
ways work with a random pure product state \ip) as the 
initial state, and p = |f/ ; )('0|. The second Renyi entropy 
S 2 (t) = — log Tr (p A (t)), like in the discussion of the main 
text for spin- 1/2 systems, is similarly related to the mea¬ 
surement of all independent operators with support in 
A: 

Tr A (p\(t)) = ^ E (Tr^W)) 2 , (13) 


Let dim("Hi) = k < 00 . Then, a basis for all Hermitian 
operators acting on a given site i is given by the trace¬ 
less, Hermitian generators T“, (a = 1, ■ ■ ■ , k 2 — 1), of the 
fundamental representation of the Lie algebra sw(fc), to¬ 
gether with the identity operator = Tf. Note that 
the space of all Hermitian operators on site i is itself a 
Hilbert space Hf of dimensionality k 2 , which should not 
be confused with the Hilbert space of quantum states 
Hi- The generators T° (a 7 ^ 0) can be chosen to obey 
the relations 


where N A is the total number of sites in A. The opera¬ 
tors O a are assumed to be Hermitian and independent, 
as defined by the Hilbert-Schmidt inner product on oper¬ 
ators: p^Tr A (( D\0’ a ) = <5 0a , 0 ^ , and they form a basis 
for all Hermitian operators on A (see Appendix A.II). 

Following the same logic as in the main text, assuming 
that each basis operator £>,4 effectively scrambles in its 
lightcone vt, we end up with: 

Tr A (p 2 A {t)) = k~^ rA+vt '> d ~ r A \ (14) 

valid until vt < ( 2 1 / d — 1 )tti, at which point Tr^p^) ~ 


rr k 2 -i 

T?T> = 6 ab h,k + I/ ^ E ^ ( 10 ) 

' ^ c—l 

where / are the antisymmetric structure constants and 
d the symmetric ones. Under the Hilbert-Schmidt inner 
product, the above operators are thus orthogonal: 

\ Tr Hi (T^T?) = \d ab Tr„, (I i|fe ) = 5 ab . (20) 
k k 

For example, for the case of k = 2, the standard basis 
operators on site i can be taken to be the familiar Pauli 
matrices: {lb, of, erf, of }. 






7 


To define a basis on the space of Hermitian operators 
on N sites T~L° = Oi'Hp, we form tensor products of T“s, 
so that a basis operator 

,T l k2 - 1 }. (21) 

There are ( k 2 ) N such basis operators in total, and they 
are orthonormal as given by the generalized Hilbert- 
Schmidt inner product: 


where 11») is the first basis vector in ~Hi for some enumer¬ 
ation of the k basis vectors (similar to | fi) in the case of 
SU( 2)) and V t the SU(k) operator (a local unitary oper¬ 
ator acting only on site i) that rotates |lj) into | fa). Let 
us call |1) = |liI 2 • • • 1 at) - Then (13) becomes 


Tr A {p\{t)) = J2^tO A U t \l) 2 : (29) 


k N 


o A 


^Tr (O'O') = ^Tr (OO') = S a ,o '■ (22) 

An example of a basis operator is O = = 

(g) T“‘ 2 <8> • • • ® T^ L ® Ia-uu-Uil which is a product 
of T“ operators on L distinct sites and a global identity 
acting on the subregion A — zi U • • • U (A is the full 
system). We say that such an operator has (physical) 
size L and has support on subregion U • • • U II- 

Since the reduced density matrix is also a Hermitian 
operator, it can be written as 

P A (t) = E c 0a (t)0 A , ( 23 ) 

Oa 

where the coefficient cq is a real number given by 

c O A (t) = ~^Tr(pO A (t)). ( 24 ) 

Squaring it gives 

P A {t)= E c o A {t)co' A {t)0 A 0' A , ( 25 ) 

0 A ,0' A 


and making use of the Hilbert-Schmidt inner product, 
(22), we get 

Tr A (p\(t)) = k N -Y, c O A (tf = E ( Tr 


o A 


o A 


(26) 


our desired result. Thus, we see that second Renyi EE 
£ 2 (t) = — logTr A (p A (t)) can be obtained by measuring, 
witin the fixed initial state p, all time evolving basis op¬ 
erators 0 A (t) that have support initially in A. 


Proving ( A ([ ) 

We want to prove 

Tr A {p\{i)) = k-^ rA+v ^ d ~ r A \ (27) 

once again for a subregion A which is a d-ball of volume 
ld,r d A = N a . 

Given a random initial product state, and following 
the main text, we can write it as 

|0i 02 ••• <Azv) = Vi | I 1 I 2 ••• in) (28) 

i 


where the unitary time evolution operator is now gener¬ 
ated by the locally rotated Hamiltonian IL^UOILK, 
which is still a local Hamiltonian having the same light- 
cone velocity as H. 

Now, we want to replace all terms (l\Uj0 A U t \l} 2 in 
the sum of |29| ) by the contribution due to a typical op¬ 
erator (l\ulO A ^y v y t \l) 2 . Following the main text, the 
average size/support of an operator in the set of basis 
operators of A is L = ( k 2 — 1 )/k 2 x N A , which follows 
from a skewed binomial distribution as there are k 2 — 1 
choices of non-trivial operators T“ acting on a site, com¬ 
pared to a choice of a single local identity Ij. In fact, an 
average operator is also typical, because the standard de¬ 
viation scales as \/N A , which is smaller than N A . Thus 
the error induced by replacing all contributions in the 
sum with that of a typical operator is small because the 
sum is dominated by the contributions from typical op¬ 
erators anyway. 

Next, we invoke the same ‘scrambling’ assumption for 
ergodic systems. Assume that a basis operator localized 
at site i scrambles in a ball of radius (vt) d ’, then, a typi¬ 
cal operator of mean size L = ( k 2 — 1 )/k 2 x N A becomes 
a complicated operator in a ball of radius ( r A + vt) d . 
Precisely, the vector of ( k 2 )' 1 A rA+vt ) coefficients co{t) 
that denotes the decomposition of a time-evolved opera¬ 
tor 0 A (t) = c o(.t)@ in terms of basis operator of the 

ball is a random vector, after time vt ~ 0(1)- 

To estimate (l\Uf OA.typ.U t \l} 2 . we have to count how 
many ‘diagonal’ operators 0 Ay diag. there are in the set of 
operators that span the ball of radius (r A + vt) d , such 
that (l|C>A,diag.|l) = 1 and (1 \O a , not diag. |1) = 0. These 
diagonal basis operators are in fact, the given by prod¬ 
ucts of T“ which live in the Cartan sub-algebra of su(k) 
(i.e. the diagonal generators), together with the identity 
operator on a site Ij. For a given site, there are k possi¬ 
ble choices only. Thus the total number of diagonal basis 
operators in the ball of radius {r A + vt) d is k ld ( rA+vt ' > . 
For the familiar case of su(2), i.e. a spin-1/2 chain, there 
is only 1 generator in the Cartan sub-algebra, of, so to¬ 
gether with the identity operator, this implies that the 
number of diagonal operators in the ball are 2 ' yd ^ rA+vt i , 
as quoted in the main text. 

Then, following the same logic as the main text, be¬ 
cause the coefficients cq in the decomposition of a time- 
evolved basis operator 0 A (t) are independent random 
variables, by assumption, the cross-terms cqCo 1 that 







arise in (1\U}OA,typ.Ut\l) 2 do not contribute, and so 


Then, 


(l\UlO A ,ty P .Ut\l ) 2 = V <£*„.(*) = 


/ -V ^^diag. 


fc'Yd(r A +vt) d 
(J i 2^'y d (r A +vt) d ‘ 

(30) 


Thus, we end up with 


TM/d(t)) ^ pvrE^l^.typ.^tl 1 ) 2 


Oa 


1 d Uld(r A +vt) d 

—(k 2 y dr *- 


kr> dr A 


(k 2 )' 1 d(r A +vt) d 


— d((r A +vt) d -r A ) 


the desired result. 


(31) 


Appendix B: Basis invariance for measurement 


Tr A (H 2 A ) = Tr A E h °- h O' A ° aO[ A 


\OaO' 


= E ho A h 0 ' A Tr A (0 A 0' A ) 
o a o' a 

= E ho A h 0 ' A k NA 5 OA o' A 

O a O' a 

= ^E^a 


o A 


k N A 


E(Tr a(HO a ))\ 


We have used O a = O a , and also the orthonormality 
of the basis set under the Hilbert-Schmidt norm in the 
above. 

Since the LHS is independent of the basis set used, 
the RHS must be as well, and the weight is given by 

k"*TrAHZ). 


(34) 


Let dim('H,) = k < oo. Given any Hermitian operator 
H , (which includes p, the density matrix), and a bipar¬ 
tition of the system into A and B , we prove that the 
‘weight’ of the operator in A, given by 

Weight = E (Tr (HOySf, (32) 

o A 

is independent of the orthonormal basis set, {O a }, that 
spans all Hermitian operators in A. The inner prod¬ 
uct on this set is given by the Hilbert-Schmidt norm: 
^Tr a { 0 \ 0 ’ a )= 8 0 a o ' a . 


II: Proof of conversation of weight of an operator 
under time evolution 

The conservation of weight of any Hermitian operator 
H under time evolution H(t) = U t HU} is a direct ap¬ 
plication of the above statement, with the subregion A 
being the full space. 

The weight of a time-evolving operator H(t) = 
ho(t)0 , given by the sum 

Em*) 2 , ( 35 ) 

o 

is conserved, since ho(t) = -p^Tr (U t HU}0) = 
-^jTr(HUjOUt), and the set {U}OU t } is an orthonor¬ 
mal basis, if {O} is, at any time t. Thus, the weight is 
constant in time and proportional to Tr (H 2 ). 


I: Proof 

The above expression can be written as 
J2o a Tr a {H a O a ) 2 , where H A = Tt b H. 

Now, H a can be decomposed as 

H A = Y J ho A 0 A . (33) 

O a 


Appendix C: Proof of replica trick 

We prove the statement 

Tr A (p n A ) = Tr n ^(nS), (36) 

where 'H® n = 1 B Ai (g> 'H.B i is an extended, compos¬ 

ite Hilbert space comprised of n copies of the original 
Hilbert space H with a bipartition into A and B. Here 
fl = p® n ? the n-replicated density matrix, where p is the 
density matrix on a single copy of the Hilbert space. S is 
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the unitary swap operator that cyclically permutes only 
states of Ai'. 


S\il 1 )A 1 1 4>2)a 2 ■ ■ ■ \i>n)A rl ® \4>)b-i,— ,B„ 

= \^n)AA^l)A 2 ■ ■ ’ It/’n-l )a„ ® \<f>)B u -,B n - (37) 


We start from the RHS, Tr w ®n(f IS). Let 

\4> ai ■ ■ ■ i’aAbi ■■■4>b n ) = It/tuUi <8 • ■ • <S> \lpa n ) A n <S> 

|</>hi)si (8» • • • <8 > |<A&„)b„> where is a set °f states 

labeled by a,; which spans Ha 4 , and similarly ) is a 


set of states labeled by 6,; which spans 'H B t ■ Then 
Tr w ®„(fi5) 

= X] ^<*1 "‘V’a„<foi ■ ■ ■ ^b n \^S\l/j ai hi 

= X! (V’ai ‘ ‘ ‘ i’a n < t>b 1 ■ ■ ■0(,Jft|V’a n ' ' ' V’an-i'K ' ' ' 06 n ) 

= X] (^ a x’--^a n \{T^Bp) m \‘lpa n ---'fpa n - 1 ) 

a 1 ,••• ,<2ri 

= ^2 (^ai\PA\^a n ){^a 2 \PA\'lpa 1 )---(i’a rl \PA\'lpa rl - 1 ) 

ai ,••• ,a n 

= i' l l’a 1 \PA\’tpa n )(i’a n \PA\^a n - 1 )---{fpa 2 \PA\fpa 1 ) 

a i ,••• ,a n 

ai 

= Tr A (p n A ), (38) 

which is the desired result. In the above, we have used 
the resolution of the identity I b a = YA a IVvX^al- 



